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Abstract 



The mischievous linear coalescence approach to hadronization of quark matter is shown to violate 
strangeness conservation in strong interactions. The simplest correct quark counting is shown to coin- 
cide with the non-linear algebraic coalescence rehadronization model, ALCOR. The non-linearity of the 
ALCOR model is shown to cancel from its simple predictions for the relative yields of (multi-)strange 
baryons. We prove, model independently, that quark degrees of freedom are liberated before hadron 
formation in 158 AGeV central Pb -I- Pb collisions at CERN SPS. 

PACS number(s): 12.38.Mh, 13.87.Fh, 24.85. +p 

1 Introduction 

It is a deep rooted desire to explain complicated experimental observations with simple and transparent 
models, understandable in laymen's terms. This very acceptable ambition inspired a recent publication 
Q, where an attempt was made to explain the relations between the multiplicities of different strange 
baryons produced in heavy ion reactions with the help of linear coalescence rehadronization model. 
This model was based on simple quark counting and elementary probability estimates. However, the 
considerations in Ref. lead to a clear violation of strangeness conservation in strong interactions: In 
order to explain the data in laymen's terms, Bialas had to assume in heavy ion reactions the number 
produced of strange quarks, s is not equal with the number of produced anti-strange quarks, s. 

We show here that this unacceptable assumption was related to the neglected requirement that all 
quarks have to hadronize, as no free quarks are observed in nature. This necessity leads to a non-linearity 
even in the simplest algebraic rehadronization model (see Ref. §). Properly taking into account the 
non-linear competition for quarks by the various hadron forming coalescence channels is a natural way to 
correct the linear coalescence treatment in the simplest possible manner. This leads to the normalized, 
non-linear a/gebraic coalescence rehadronization model, ALCOR |^], which automatically takes care of 
the conservation of strangeness in strong interactions as well. 

2 Linear coalescence in rehadronization? 

In the linear coalescence rehadronization model Bialas assumed, that the number of produced particles 
is proportional to the product of the numbers of constituent particles within the reaction volume. In 
the following we denote the number of particles by the particle symbols, in particular, q stands for the 
number of light (up and down) quarks. Linear coalescence yields the following relations: 



where A|E = A -I- E stands for the total number of strange baryons that contain a single strange quark. 
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Similarly, for the anti-baryons one obtains 
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also assumed, that the coefficients of proportionality for particles and their antipar- 
etc. As the coefficients a^, an describe effectively the (particle anti-particle 
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Furthermore, Ref. 
tides are equal: a 

symmetric) likelihood that a multi-quark bound state is formed once the constituents are given, this 
seems to be a very reasonable model at first sight. 

This linear model is attractive not only because of its simple formulation but also because of its 
simple predictions, namely that the unknown a coefficients cancel from the anti-baryon to baryon ratios. 
This leads to the following equations: 
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In order to explain the value of 



= 0.383 ±0.081, 
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as measured in Pb -|-Pb reactions at CERN SPS energies, Bialas had to assume s 7^ s as a consequence of 
eq. ( p^ of his linear model, which is a clear violation of strangeness conservation in strong interactions. 
What is the origin of this contradiction? 

The linear coalescence model, as defined in eqs. (3) and (4) is a good approximation only tf the 
composite particles use up a small fraction of the constituents. This is the case e.g. in the coalescence 
treatment of the deuteron (j>n) or triton (pnn) formation from the gas of protons and neutrons. In that 
situation most particles will remain in nucleon (p or n) state. Thus the competition is negligible among 
the tritons and the deuterons for the building block nucleons in that case. 

In the case of hadronization of a quark matter all of the constituent particles (the quarks) have to be 
placed into composite particles, namely into colorless hadrons. No free quarks are observed in Nature. 
This is the very essence of color confinement. Thus, during the hadronization process, a fixed number of 
strange quarks must be distributed among hyperons and anti-K mesons. The anti-strange quarks have 
to be distributed among anti-hyperons and K mesons. [ Note that K meson is a (gs) bound state, while 
the anti-K meson is a (gs) bound state.] However, in heavy ion collisions we have incoming nuclei in 
the initial state, that contain nucleons, formed by constituent quarks, only. As quarks are produced in 
quark antiquark pairs, the number of quarks in the final state is larger than the number of anti-quarks. 
Hence we must have more K mesons than anti-K mesons. 

The processes creating different hadrons are not independent, they compete with each other, contrary 
to the basic assumption in the linear coalescence model. 

The redistribution can be counted for by introducing normalization factors, bs,h-^,hq,hq. These bi 
normalization factors are not free parameters! They are determined from the requirement that all quarks 
must be recombined into hadrons during the hadronization process, as described first in the ALCOR 
model, Ref. [i. 
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3 The ALCOR approach 
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The non-linear coalescence equations for the formation of quark-antiquark clusters during rehadronization 
read as: 

(14) 
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The meson yields are determined similarly, 

n'' = C^b.b^qq, (22) 

K = CKb.b^qs, (23) 

K = Cj^b^bsqs, (24) 

rj = Cr,bsb^ss. (25) 

The normalization coefficients bq, bs, bq and b^ are determined uniquely by the requirement, that the 
number of the constituent quarks do not change during the hadronization — which is the basic assumption 
for all quark counting methods: 

s = 3f7 + 2H + A|E + ^+ 77, (26) 
s = 3n + 2H + A|S + A' + T), (27) 
q = 3p + 2A|E + H + i^ + 71^ (28) 
q = 3p + 2A|E + H + ¥+7r'* . (29) 

Here n'^ is the number of directly produced [qq) states. ( Note that other particle symbols also stand for 
the number of directly produced particles of a given type, and care must be taken when comparing the 
ALCOR predictions with data, especially regarding the corrections of particle numbers for the feed-down 
of resonance decay contributions.) 

Similar, non-linear coalescence equations with constraints of using up all the coalescing particles were 
proposed in a simpler form first in Ref. in the so called combinatoric break-up model, assuming a 
constant coalescence coefficient for all the baryons and another one for all the mesons. Later the ALCOR 
model was formulated as presented above in a simple form. ALCOR equations were solved numerically 
in Refs. |^ ^, ^ , using, as input values to the coalescence equations, the number of quarks as predicted 
by Monte Carlo event generators. The ALCOR equations were derived from a set of rate equations in 
Ref. [^, in the sudden ap pro ximation. 

Substituting eqs.(p^ - |25|) into the constraints given by eqs. ( p6| - p^ one obtains a set of non- linear 
equations for the bi normalization constants. However, one can predict some relations even without 
solving these set of nonlinear equations. It turns out that the particle-antiparticle ratios can be expressed 
in terms of the effective number of quarks even in the non-linear ALCOR model, similarly as it was done 
by Bialas in the linear coalescence model. 

Let us denote by upper case Q, S the reduced, effective number of up-|-down and strange quarks, in 
contrast to the lower case q and s, that stand for the total number of (up+down) and strange quarks 
before the hadronization: 

Q = bqq, (30) 
Q = (31) 
S = b,s, (32) 
5 = (33) 
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Furthermore, one could make the usual assumption that the C coalescence coefficients for baryons are 
equal to that of the corresponding antibaryons, 



Ca = C^, 
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With these notations, the following relations are obtained in the non-linear ALCOR model: 
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These equations are formally similar to the equations of the linear treatment displayed in eqs. (p[]l2[), 
however, they are given in terms of the effective number of quarks, Q = bqq, that are complicated, 
non-linear functions of the number of quarks q, s available before the hadronization. Thus one easily gets 
the following interesting relations: 
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These are the relations among the ratios of the observable number of particles that should be satis- 
fied if the particle production in some reaction proceeds via algebraic recombination of the independent 
valence quarks. This way a direct relation is obtained in a self-consistent manner, that connects the 
presence/absence of independent, initial quarks with the observable yields of hadrons. This relation is 
to a large extent model-independent (i.e. independent of the initial quark content, independent from the 
values of the coalescence coefficients C and independent from the values of the non-linear renormaliza- 
tion factors bi). Only two physical assumptions enter these model-independent eqs. ( ]4^ - |45[ ): that the 
rehadronization process is sudden and that the valence quarks are available in unbound states before the 
hadronization. 

One may evaluate the following measurable numbers: 
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If hadronization proceeds via a sudden recombination of quarks, then all these d numbers should be 
equal, d = (Ia = ds = dn. 

The experimental values for the anti-baryon to baryon number ratios are determined in Refs. ^, ^, ^ 
as follows: 



N 
A|E 



K+ 



0.070 ±0.010, (50) 

0.133 ±0.007, (51) 

0.249 ±0.019, (52) 

0.383 ± 0.081, (53) 

d = 1.80 ±0.2 (54) 



K- 

which yields the following values for the d factors of strange baryons : 

dA = 1.9 ±0.3, (55) 
ds = 1.89 ±0.15, (56) 
dn = 1.76 ±0.15. (57) 

These numbers show a very good agreement with d, a statistically acceptable, good agreement between 
the ALCOR model and the experimental data in case of Pb + Pb reaction at CERN SPS. This indicates, 
regardless of the details of the confinement mechanism, that hadron production in Pb +Pb reaction at 
CERN SPS proceeds via a sudden and complete coalescence of constituent quarks to hadrons. 

Taking into account the conservation of strangeness in strong interactions, which demands 1 — s, 
from the fi/Sl ratio we arrive at 

^ = ^ = 0.75 ± 0.06. (58) 

O Us 

Thus we obtained agreement with the experimental data without assuming s s. The deviation of 
Q/Q, from unity is caused by the difference in the normalization factors 6s and 65-, which can easily be 
understood: there are more quarks then antiquarks in the initial system, and thus more s quarks are 
used up in the production then s quarks in the K~ meson creation. Thus less s remains for the SI 
production then s quarks for the Q production. 

In the time-dependent solution of confining rate equations for the hadronization jll|, conservation 
laws are ensured by the structure of the coupled system of differential equations. The ALCOR model, 
as presented above, can be reobtained from a set of non-linear rate equations in the sudden approx- 
imation [^, in the limit when the time of hadronization is very short. Indeed, a detailed analysis of 
single particle spectra and two-particle correlation data indicates a short duration. At ~ 1.5 fm/c, for 
the production of final state hadrons at CERN SPS Pb + Pb reactions jlo|. This result excludes a 
long-hved, evaporative mixture of quark-gluon plasma and hadronic phase, that would produce pions 
over an order of magnitude larger period. Thus the experimental data prefer a sudden production of 
hadrons, a process that may happen to be out of local thermal and chemical equilibrium. 

Finally, let us point out, that the strangeness conservation leads to the following relation: 



^ n H „ S K ^ AIE 

3— — Tj-±2— — ±— — ±1 = ' 



H A|S A|E A|E AjE 



3^^ + 2^ + JL + l 
H A|E A|E A|E 



(59) 
(60) 



With the SPS Pb+Pb data |^, ^ the left hand side of this equation is 2.57, while the right hand side is 
2.66. Thus the strangeness conservation equation is also well fulfilled for this case. The validity of eq. ( |6o| ) 
is an absolute measure of strangeness conservation in strong interactions. If the experimental values do 
not satisfy eq. (|60|), then the detectors miss important parts of the particle momentum distributions, or 
they have a serious systematic error. 

Further observations of Ref. namely dA = ds for the SPS ^^S +^'^ S reactions, but cJa 7^ ds for 
the SPS p±Pb reaction remain unchanged with our analysis. 



5 



4 Conclusions 



We clarified the reason why the hnear coalescence model should not be applied to the hadronization of 
the quark matter: it leads to a violation of strangeness conservation. 

We have shown, however, that this shortcoming can be corrected if a a nonlinear coalescence model 
is introduced, that takes into account the conservation of strangeness in strong interactions and the 
fact that no free quarks are observed in Nature. In particular, we have shown that the solution of 
the ALCOR hadronization model yields simple, parameter independent relations in a consistent manner. 
These relations were surprisingly simple and similar to those obtained from the linear coalescence picture. 
Data on strange particle production in Pb+Pb reactions at CERN SPS satisfy the coalescence model 
predictions, while they are known to be not satisfied by data from p + Pb reactions at CERN SPS. 
Based on considerations of the production ratios of strange hadrons, we have obtained an elementary, 
self-consistent and model independent proof, that quark degrees of freedom are liberated in Pb + Pb 
reactions at CERN SPS before the onset of hadron formation. 
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